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Intensity-modulated particle therapy (IMPT) with carbon ions is comparatively susceptible to various 
uncertainties caused by breathing motion, including range, setup, and target positioning uncertainties. To deter- 
mine relative biological effectiveness-weighted dose (RWD) distributions that are resilient to these uncertain- 
ties, the reference phase-based four-dimensional (4D) robust optimization (RP-4DRO) and each phase-based 4D 
robust optimization (EP-4DRO) method in carbon-ion IMPT treatment planning were evaluated and compared. 
Based on RWD distributions, 4DRO methods were compared with 4D conventional optimization using planning 
target volume (PTV) margins (PTV-based optimization) to assess the effectiveness of the robust optimization 
methods. Carbon-ion IMPT treatment planning was conducted in a cohort of five lung cancer patients. The 
results indicated that the EP-4DRO method provided better robustness (P = 0.080) and improved plan quality 
(P = 0.225) for the clinical target volume (CTV) in the individual respiratory phase when compared with 
the PT V-based optimization. Compared with the PTV-based optimization, the RP-4DRO method ensured the 
robustness (P = 0.022) of the dose distributions in the reference breathing phase, albeit with a slight sacrifice 
of the target coverage (P = 0.450). Both 4DRO methods successfully maintained the doses delivered to the 
organs at risk (OARs) below tolerable levels, which were lower than the doses in the PT V-based optimization 
(P < 0.05). Furthermore, the RP-4DRO method exhibited significantly superior performance when compared 
with the EP-4DRO method in enhancing overall OAR sparing in either the individual respiratory phase or refer- 
ence respiratory phase (P < 0.05). In general, both 4DRO methods outperformed the PT V-based optimization 


in terms of OAR sparing and robustness. 
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I. INTRODUCTION 


Intensity-modulated particle therapy (MPT) with carbon 
ions can deliver higher doses to the target volume while 
significantly sparing the adjacent organs at risk (OARs) 
[1, 2]. The sharp dose falloff behind the pristine Bragg 
peak of carbon-ion beams makes the dose distributions highly 
sensitive to uncertainties, which greatly diminishes the effec- 
tiveness of the IMPT [3]. Three of the most pertinent uncer- 
tainties in carbon-ion dose distributions include (1) the range 
uncertainty introduced by the uncertainty in computed to- 
mography (CT) numbers and their conversion to relative stop- 
ping power [4]. (2) Setup uncertainty is related to a lack of 
reproducibility in patient positioning [5]. (3) uncertainties in- 
troduced by respiratory motion, especially for patients with 
lung cancers [6]. The comparison of PTV-based optimization 
with robust optimization is a widely used approach for eval- 
uating the robustness of treatment plans. Although the PTV 
concept has been employed in some studies of carbon-ion ra- 
diotherapy [7], it is insufficient for carbon-ion radiotherapy. 
Dose distributions are influenced by various factors, not only 
at the edges of the target volume but also inside the target 
volume [8]. 
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To mitigate the impact of various uncertainties in carbon- 
ion IMPT, various three-dimensional robust optimization 
(3DRO) methods have been developed, including probabilis- 
tic optimization [9, 10], voxel-wise worst-case robust opti- 
mization [11], and worst-case scenario robust optimization 
[12]. These methods aim to incorporate uncertainties directly 
into optimization. Probabilistic methods optimize treatment 
plans based on a large number of dose distributions produced 
by a random sampling setup and range uncertainty scenar- 
ios with an assumed probability [3]. Voxel-wise worst-case 
robust optimization considers the minimum and maximum 
doses from all uncertainty scenarios for each voxel inside the 
target volume and maximum doses for each voxel inside the 
normal tissues. The worst-case dose distributions include the 
minimum and maximum doses. During the iterations, only 
one objective function based on the worst-case dose distri- 
butions was optimized. The worst-case scenario robust opti- 
mization method evaluates the objective function for all un- 
certainty scenarios and chooses the worst objective function 
score in the iteration process. Recent studies have shown that 
3DRO methods can minimize the variance in dose distribu- 
tions under different uncertainty scenarios and improve the 
robustness of treatment planning [13-17]. 


Target motion and motion-induced range changes play im- 
portant roles in particle radiotherapy. The effectiveness of 
3DRO methods in mitigating the effect of respiratory motion 
on lung cancer therapy is limited [18]. To explicitly account 
for respiratory motion in IMPT, different strategies have been 


proposed, mainly including respiratory gating and breath- 
holding. Unfortunately, these methods have limitations in 
terms of treatment time and patient tolerance requirements. 
Additionally, online motion tracking [19] is indispensable for 
moving targets in the treatment process to reduce the im- 
pact of breathing motion. However, this remains technically 
challenging. Another approach proposed by Graeff et al. [20] 
suggested that 4D optimization based on all motion states 
is valid for carbon-ion therapy. Liu et al. [21] introduced 
a reference phase-based 4D robust optimization (RP-4DRO) 
method, which optimizes cumulative 4D dose distributions. 
Ge et al. [3] proposed a method termed as phase-based 4D 
robust optimization (EP-4DRO), which optimizes dose dis- 
tributions based on individual respiratory phase. Wolf et al. 
[7] proposed a robust nonlinear RBE-weighted optimization 
method to expand the carbon-ion IMPT, thus exploring the 
potential for improving plan robustness and sparing critical 
organs. These studies confirm that 4DRO methods can im- 
prove plan robustness and offer better control over uncer- 
tainties during treatment. However, there has been a persis- 
tent lack of comparative analyses among the different 4DRO 
methods in carbon-ion IMPT treatment planning, hindering 
a comprehensive understanding of the distinct characteristics 
of each method. 


The relative biological effectiveness (RBE), which is de- 
fined as the ratio of a photon dose and an ion dose leading to 
the same biological effect under identical conditions [22], is 
widely used in particle therapy. A constant RBE value of 1.1 
was clinically accepted for protons. 


However, the RBEs of carbon ions must be described by 
mathematical models rather than by a single parameter be- 
cause of their RBE values’ sophisticated dependence on phys- 
ical and biological factors. It is necessary to plan carbon-ion 
treatment with RBE-weighted doses (RWD), defined as the 
product of physically absorbed doses and RBE. Some mod- 
els have been proposed, such as the mixed beam model, the 
local effect model (LEM), and the microdosimetric kinetics 
model (MKM) [22]. These models facilitate carbon-ion treat- 
ment planning using RWD. 


An appropriate robust optimization method is crucial for 
addressing the dose perturbations arising from various uncer- 
tainties. To the best of our knowledge, no previous studies 
have compared the RP-4DRO and EP-4DRO methods. It 
is acknowledged that the RP-4DRO method does not explic- 
itly optimize the dose distributions of the individual breath- 
ing phases. It is necessary to pay attention to the dose dis- 
tributions in individual breathing phases for clinical imple- 
mentation. To evaluate the effectiveness of 4DRO methods 
in carbon-ion IMPT, a comparative and analytical investiga- 
tion of 4DRO methods is indispensable. The innovation is to 
compare and evaluate two theoretically different robust opti- 
mization methods in terms of RWD distributions in carbon- 
ion IMPT treatment planning. 


Il. MATERIALS AND METHODS 


A. Treatment planning 


MatRad [23], an open-source treatment planning system 
(TPS), was developed for educational and research purposes. 
This also supports carbon-ion IMPT treatment planning. The 
energies of the carbon ions available in matRad are in the 
range of 115.23 MeV /u to 398.84 MeV /u corresponding to 
a penetration range of 32.68 mm to 294.25 mm in water. Fur- 
thermore, alpha-beta ratios of 10 and 2 were selected for the 
target volume and normal tissue, respectively. The dose grid 
was set to 3mm Xx 3mm x 3 mm. The phosphate buffered 
saline (PBS) method with spot scanning was used. Further- 
more, PBS can realize conformal dose distributions for irreg- 
ularly shaped tumors and improve the beam utilization [24]. 


Gaussian-shaped carbon-ion beam spot sizes range from 
11.80 mm to 8.01 mm (full width at half maximum, FWHM) 
for energies in the range of 115.23 MeV /u to 398.84 MeV /u. 
A spot spacing of 3 mm in the lateral and longitudinal di- 
rections was employed. Carbon-ion IMPT treatment plan- 
ning was performed in five patients with lung tumor le- 
sions. For each patient, 4DCTs were acquired with 10 motion 
states. To reduce the consumption of computer memory, both 
4DRO methods were developed using three significantly dif- 
ferent breathing phases: 0% expiration phase, 80% expiration 
phase, and 40% inspiration phase. The 0% expiration phase 
is defined as the reference breathing phase. 


For the PTV-based optimization, the PTV was generated 
by extending the union of the CTVs in all three phases with a 
5-mm safety margin. The magnitude of the margins was de- 
rived from the original photon treatment planning [25]. The 
dose distributions were calculated for the three respiration 
phases without accounting for setup and range uncertainties. 
Both 4DRO methods were performed on the CTV for each 
phase without incorporating any margins. The setup uncer- 
tainties were assumed as +5 mm in +z, +y, and +z di- 
rections. The range of the uncertainty parameters was set 
to +3.5% when the water-equivalent depth was computed. 
Each breathing phase had eight uncertain dose scenarios and 
one nominal dose scenario. The tumor size, location, mag- 
nitude of motion, and beam angle of the five patients un- 
der investigation are listed in Table 1. The CTV centroids 
in the right-left (RL), anterior—posterior (AP), and supe- 
rior—inferior (SI) directions for the three breathing phases 
were recorded. The maximum centroid is used to deter- 
mine the magnitude of the motion. The patient datasets in- 
cluded 4DCT images with CTV and OARs delineated by 
an experienced physician following the delineation guide- 
lines provided in [26]. The prescription for all patients was 
60 Gy(RBE)/15 fractions with threshold doses of 45, 20, 
and 30 Gy(RBE) for the spinal cord, lungs, and heart, re- 
spectively. The LEM was used to calculate RBE values. 


TABLE 1. Tumor locations, volumes, and motion magnitudes of the 
patients under investigation. 


Tumor Motion Beam 
7 Tumor 
Patient ID À volume mag. angle 
location = 
(cc) (mm) C) 
1 Right hilar 47.23 20.0 220 / 140 
2 Right hilar 56.99 19.0 220/ 140 
3 Left hilar 40.86 12.0 210/ 140 
4 Right lower 1.33 2.0 330 / 45 
5 Right lower 64.31 5.0 330 / 45 
B. 4D robust optimization 
Liu et al. [21] proposed a robust version of the 4D- 


optimization method and demonstrated that this approach can 
enhance the robustness of the cumulative 4D dose distribu- 


tion. The physical absorbed dose distributions (dj, jw?) in 
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different setups and range uncertainty scenarios were calcu- 
lated for the breathing phases of 4DCT datasets, where wj 
represents the intensity weight of the beamlet j, and dj, j 
represents the influence matrix describing the distribution of 
the beamlet j in the breathing phase k to voxel i, in the un- 
certainty scenarios. We define s as each uncertainty scenario 
considered during treatment planning optimization and define 
S as the set of all scenarios. During the optimization process, 
a non-negative vector w? is used to transform the constrained 
optimization problem (w > 0) into a non-constrained prob- 
lem. The cumulative physically absorbed dose distributions 
were deformed to the reference breathing phase by using the 
registration algorithm in matRad. The cumulative RWD dis- 
tributions D? are computed using the cumulative physically 
absorbed dose distributions. Specifically, D? are employed in 
a standard quadratic objective function to design a carbon-ion 
treatment plan. In the objective function, the worst-case sce- 
nario robust optimization mentioned above was conducted as 
follows: 
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where Dpre and D, denote the prescription dose in the tar- 
get volume and dose constraints on the relevant tissues, re- 
spectively. Furthermore, N and M denote the numbers of 
target volumes and OARs, respectively. Parameters qn maz» 
dn,min» Nd qm,max are the penalty factors for the target vol- 
umes and OARs. Furthermore, qn,maz = 1.0, qn,min = 1.0, 
and qm,max = 1.0 were used. The step function H is math- 
ematically defined as follows: H = 1.0 when D? > Do and 
H = 0.0 otherwise. The ultimate objective of the optimiza- 
tion is to minimize the objective function f(w) to realize the 
optimal dose distributions over the relevant voxels. 
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The RP-4DRO method did not explicitly address the op- 
timization of dose distributions in the individual breathing 
phases ; it remained unclear whether robustness in individ- 
ual respiratory phases can be ensured. Ge et al. [3] proposed 
the EP-4DRO method for solving the uncertainty problem in 
proton treatment planning. The RWD distributions D?, in dif- 
ferent setups and range uncertainty scenarios were calculated 
for all the breathing phases of 4DCT. The worst dose distri- 
butions for each breathing phase were calculated by selecting 
the worst objective score. Subsequently, the worst dose dis- 
tributions for all the phases were simultaneously optimized. 
The objective function is expressed as follows: 
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It should be noted that all the parameter settings were the 
same as those in the RP-4DRO method. During the optimiza- 
tion iteration, the objective function values for each dose sce- 
nario were calculated and the worst scenario was selected. 


C. Plan evaluation 


To maintain consistency, the resulting dose distributions for 
the volume of interests (VOIs) were evaluated and compared 
using the same parameters. Nine dose distributions for each 
breathing phase (6 setup uncertainty scenarios, 2 range un- 


certainty scenarios, and 1 nominal scenario) were calculated. 
The bands of dose-volume histograms (DVHs) derived from 
the 27 dose distributions (nine scenarios for three breathing 
phases) were analyzed to compare the robustness of the treat- 
ment plans. For PTV-based optimization, the dose distribu- 
tions were recalculated for each breathing phase using the op- 
timized set of spot intensities for each of the nine uncertainty 
scenarios. A narrower band in the DVH indicates greater ro- 
bustness, implying that the dose distribution changes are rel- 
atively minimal under the influence of uncertainties. Evalu- 
ating DVHs with nine accumulated dose distributions is es- 
sential [3]. Thus, DVHs derived from the 9 accumulated dose 
distributions were also analyzed. Dose accumulation was as- 
sessed during the three breathing phases. 

Dose distribution conformity index (CI) and target hetero- 
geneity index (HI) were computed using the following equa- 
tions: 
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Specifically, TV35% denotes the target volume receiving 
a 95% prescription dose. Furthermore, V35% represents the 
total volume receiving at least 95% of the prescribed dose; 
CTV is the target volume for the three methods. The CI is a 
measure of how well the target dose distributions conform to 
the target volume, with a value closer to 1.0 indicating greater 
conformity. HI reflects the homogeneity of the target dose 
distributions, with a value closer to 0.0 indicating greater ho- 
mogeneity. To compare the target coverage, the index Dos, 
is used to represent the dose received by the 95% target vol- 
ume. A value closer to the prescribed dose indicates better 
TV dose coverage of the target volume. 

Normal tissue sparing was evaluated by comparing differ- 
ent treatment plans optimized using different methods. Evalu- 
ation indexes, such as lung Dmean, lung V 20, lung V5, heart 
Dmean, heart V30, and spinal cord Daz were analyzed and 
compared among the aforementioned methods. Given that 
the dose distributions were non-Gaussian profiles, a statisti- 
cal analysis was conducted using the Wilcoxon signed-rank 
test [27] to compare the PTV-based optimization and 4DRO 
methods. A statistically significant difference was considered 
when the P value was less than 0.05. 


Il. RESULTS 
Five lung cancer patients were evaluated for potential dif- 


ferences. All methods focused on RWD. This section presents 
key findings pertaining to these patients. 


A. Exemplary patient 


Fig. 1 shows the DVH bands for the CTV and OARs de- 
rived from the 27 dose distributions for Patient 1 as an exam- 
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Fig. 1. DVHs of Patient 1 for the different treatment plans using the 
PTV-based optimization (a), RP-4DRO (b), and EP-4DRO (c) meth- 
ods. The axial view of the accumulated 4D dose distribution overlaid 
on the reference phase is displayed on panels (d—f), corresponding 
to the different plans. 


ple. The results of the PTV-based optimization, RP-4DRO 
and EP-4DRO methods are shown in Fig. la, 1b, and Ic 
respectively. The dose distributions without considering the 
uncertainties for the reference phase ( nominal scenario) are 
highlighted by black solid lines. For Patient 1, the bandwidths 
at Dgsy, for the PTV-based optimization, RP-4DRO, and EP- 
4DRO methods were 8.92, 8.77, and 5.29 Gy(RBE), respec- 
tively. Correspondingly, the bandwidths at D5% were 2.11, 
2.11, and 1.68 Gy(RBE). The planning objective of CTV cov- 
erage was set to Dg5%/Dpre > 95% in the nominal sce- 
nario. The CTV coverages in the nominal scenario for the 
PTV-based optimization, RP-4DRO, and EP-4DRO methods 
were 94.86%, 93.51%, and 96.73%, respectively. 

For the DVHs of the normal structures, both 4DRO meth- 
ods realized narrower bands when compared to the PTV- 
based optimization in lung, heart, and spinal cord spar- 
ing. A nominal view of the accumulated 4D dose distri- 
butions in the reference phase for Patient 1 is presented in 
Fig. ld, le, and 1f. The dose at the edges of the CTV for 
both 4DRO methods was reduced when compared to the PTV- 
based optimization. Notably, the RP-4DRO method (Fig. le) 
showed hot spots in the medial area of the CTV when com- 
pared with the EP-4DRO method (Fig. 1f). 


B. Entire patient cohort 


Fig. 2 shows the target coverage indices for the nominal 
scenario. 

Compared with the PTV-based optimization, the target 
coverage indices realized by the RP-4DRO method were rel- 
atively low, except for Patient 3 (P = 0.383). The target 
coverage indices were also below constraint (Dos, / Dore = 
95%). Conversely, the EP-4DRO method exhibited supe- 
rior performance in terms of target coverage for all patients 
(P = 0.225). 


100 m PTV-based 
99 = RP-4DRO 
98 m EP-4DRO 


a 
a 


Target Coverage (%) 
ooo Oo OO 
e Oo 
L] 
a 


wo 
a 


3 
Patient ID 


Fig. 2. Target coverage indices for the nominal scenario. 


Fig. 3 shows the robustness evaluation results and plan 
qualities with three optimization methods from five patients. 
The bandwidths at 95% and 5% CTV volumes for all three 
methods are shown in Fig. 3a and 3b, respectively. A lower 
value indicated a higher level of robustness. Based on the 
experimental findings, the dose distributions within the CTV 
derived from EP-4DRO exhibited narrower bandwidths than 
those achieved through PTV-based optimization, except for 
patient 3. Similarly, the RP-4DRO method yielded narrower 
bandwidths in the dose distributions when compared to the 
PTV-based optimization, with the exception of patients 2 and 
3: 

Fig. 3c and 3d present the values of HI and CI for all 27 
scenarios across the five patients. Error bars have been added 
to the figure to indicate one standard deviation of the results. 
For the EP-4DRO method, the HI values were generally lower 
than those for the RP-4DRO method and PTV-based opti- 
mization (P = 0.345), except for Patient 4. Lower values 
indicate better dose homogeneity at the target boundaries. 
Conversely, the RP-4DRO method yielded relatively higher 
HI values, indicating worse dose homogeneity in the individ- 
ual respiration phases (P = 0.107). As shown in Fig. 3d, 
the CI values obtained using the EP-4DRO method are higher 
compared to the PTV-based optimization (P = 0.160) and 
RP-4DRO method. The CI values produced by the EP-4DRO 
method were higher than those of the PTV-based optimization 
except for Patient 4 (P = 0.005). 

As shown in Fig. 4, the dose-volume indices for the 
OARs are also calculated. A comparative analysis revealed 
that the RP-4DRO method delivered lower doses to the lung 
(deft/right) lung Dean, P < 0.05; left lung V20, P = 
0.068; (left/right) lung V5, P < 0.05; and right lung V 20, 
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Fig. 3. Bandwidth at 95% CTV volume (a) and bandwidth at 5% 
CTV volume (b) for the PTV-based optimization, RP-4DRO, and 
EP-4DRO approaches. The error-bar of heterogeneity (c) and con- 
formity indices (d), respectively. 
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Fig. 4. Dose distributions in lung for the PT V-based optimization, 
RP-4DRO, and EP-4DRO methods. 


P < 0.05). As shown in Fig. 5, both 4DRO methods yielded 
lower indices for the heart and spinal cord when compared to 
the PTV-based optimization (P < 0.05), indicating the heart 
and spinal cord received less doses. Notably, for the five pa- 


tients, the RP-4DRO method consistently exhibited the lowest 
doses to the lungs and heart across the individual breathing 
phases. However, the RP-4DRO method did not demonstrate 
superiority in maintaining lower doses delivered to the spinal 
cord when compared with the EP-4DRO method. The P val- 
ues indicate that the differences in the OARs indices between 
the two 4DRO methods were statistically significant. 
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Fig. 5. Dose distributions in heart and spinal cord for the PTV-based 
optimization, RP-4DRO, and EP-4DRO methods. 


The aforementioned results illustrate 27 dose distributions 
for the three respiratory phases. However, for further analysis 
in a clinical implementation, it is necessary to compare and 
calculate the accumulated dose distributions over the three 
respiratory phases. 

The bandwidth at 95% CTV volume and target coverage 
of the cumulative distributions are shown in Fig. 6. The 
RP-4DRO method showed a smaller bandwidth at 95% CTV 
volume when compared to the PTV-based optimization (P = 
0.022). Conversely, the EP-4DRO method does not guarantee 
small bandwidths for Patient 2 and Patient 5 when compared 
with the PTV-based optimization (P = 0.272). The target 
coverages generated by the RP-4DRO method (P = 0.45) 
and EP-4DRO method (P = 0.066) were lower than that of 
the PTV-based optimization. The dose distributions in OARs 
are illustrated in Fig. 7. Compared with the PTV-based op- 
timization, both 4DRO methods achieved lower doses in the 
lung, heart, and spinal cord (P < 0.05). Notably, the EP- 
4DRO method resulted in the lowest doses in the lungs and 


heart. In contrast, the EP-4DRO method was not superior in 
maintaining the lowest cumulative dose distribution delivered 
to the spinal cord. This implied that the accumulated dose 
distributions generated by both 4DRO methods are more con- 
servative. 
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Fig. 6. Accumulated results of robustness and target coverage for 
the PTV-based optimization, RP-4DRO, and EP-4DRO methods. 
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Fig. 7. Accumulated dose distributions in lung, heart, and spinal 
cord for the PTV-based optimization, RP-4DRO, and EP-4DRO 
methods. 


IV. DISCUSSION 


Both 4DRO methods have been proposed to improve the 
robustness of treatment plans for carbon-ion IMPT. Although 
many studies have investigated the responses of patients with 
diverse cancers to various optimization methods, there have 
been no direct comparisons between the RP-4DRO and EP- 
4DRO methods. Consequently, this study aimed to systemat- 
ically compare and evaluate both 4DRO methods in carbon- 
ion IMPT treatment planning. 

The EP-4DRO method can ensure plan robustness in indi- 
vidual respiratory phases, as shown in Fig. 3a and 3b. Con- 
versely, the RP-4DRO method ensures the robustness of the 
cumulative dose distributions, specifically in the reference 


respiratory phase, as illustrated in Fig. 6a. This is because 
the RP-4DRO method does not explicitly optimize the dose 
distributions of individual breathing phase. 


Regarding target heterogeneity (Fig. 3c and 3d), the EP- 
4DRO method shows superior dose homogeneity inside the 
target boundaries when compared to the PT V-based optimiza- 
tion and RP-4DRO methods. Both 4DRO methods outper- 
formed in the dose distribution conformity when compared 
to the PTV-based optimization except for Patient 4. This dis- 
crepancy can be attributed to the inherent challenge of en- 
suring dose-distribution conformity for smaller target move- 
ments during the optimization process. 


It is vital to guarantee the target coverage for carbon-ion 
IMPT treatment planning. Regarding the target coverage 
indices in the individual respiratory phases, the EP-4DRO 
method ensured the target coverage. Conversely, the RP- 
4DRO method sacrificed the target coverage, except for pa- 
tient 3. As shown in Fig. 6b, both 4DRO methods sacri- 
ficed the target coverage in the reference respiratory phase. 
Notably, the target coverage remained unsatisfactory for the 
uncertainty scenarios, as shown in Fig. 1. Furthermore, the 
target coverage achieved by the EP-4DRO method in the in- 
dividual respiratory phases differed from that observed in the 
reference phase. This might be due to the accuracy of data 
processing such as deformable image registration. 


Both 4DRO methods demonstrate superior OARs sparing 
when compared with the PTV-based optimization, as indi- 
cated in Fig. 4, 5, and 7. An important finding demonstrated 
that the RP-4DRO method exhibited superior performance to 
the EP-4DRO method in protecting the lungs and heart in 
both the individual respiratory phase and reference respira- 
tory phase. However, the RP-4DRO method did not exhibit 
any advantage over the EP-4DRO method in terms of spinal- 
cord sparing. The P-values indicate that the differences in 
OARs indices between the two methods were statistically sig- 
nificant. 


Generally Both 4DRO methods yielded different results 
because of the different dose distributions selected during the 
optimization process. Ge et al. [3] demonstrated that the EP- 
4DRO method maintained robustness even if the patient did 
not breathe consistently, fraction-by-fraction, over the treat- 
ment course. Our study suggests that carbon-ion IMPT treat- 
ment planning using the EP-4DRO method is superior to the 
RP-4DRO method in terms of robustness and planning qual- 
ity. while being flexible. The RP-4DRO method exhibited 
superior OARs protection. Thus, if we focus on OAR spar- 
ing, then the RP-4DRO method should be used. Conversely, 
if robustness is a priority, then the EP-4DRP method should 
be chosen. Furthermore, the results highlight that carbon-ion 
IMPT treatment planning with both 4DRO methods lead to 
worse target coverage for a small target volume, as observed 
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V. CONCLUSIONS 
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